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.  A  three-dimensional  parabolic  Navier-Stokes  numerical  solution  algorithm  has  been  analyzed  for 
prediction  of  the  nearfield  flow  development  of  a  VSTOL  Jet  in  subeonic  croaeflow.  The  esssn- 
tie!  aspects  of  algorithm' definition,  with  regards  to  initial  and  boundary  condition  pacifications, 
has  bean  summarized.  A  momentum  conserving  pressure  gradient  computation  has  baen  devel¬ 
oped  to  complete  the  problem  definition  and  facilitate  problem  initialization.  A  sequence  of 
overlapping  interaction  solutions  has  been  evaluated  for  prediction  of  a  circular  Jet  at  X  ■  8.0. 
The  results  of  the  3DPNS  predictions  have  compared  qualitatively  with  the  sparse  available  ex¬ 
perimental  data.  It  is  crucial  that  quality  data  be  acquired  to  permit  quantitative  assessment  of  the 
results  of  this  analysis  6rocadure.  ^ _ _ _ 
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INTRODUCTION 

Analysis  and  experimental  evaluation  of  the  aerodynamics  of  a  subsonic 
VSTOL  .iet  in  cross-flow  has  been  a  major  focus  of  a  segment  of  the  technical 
community,  cf.,  references  [1-7].  Figure  1  is  a  schematic  illustration  of 
the  basic  geometry.  The  fundamental  differences  in  the  cross-flow  velocity 
field  immediately  adjacent  to  plane  of  injection  of  the  VSTOL  jet,  Z  =  0  in 
Figure  1,  and  a  solid  cylinder  of  the  same  diameter,  is  evident  in  the  oil- 
flow  streakline  data  of  McMahan  and  Mosher  [8],  see  Figure  2.  Figure  3,  a 
close-up  view  of  similar  data  obtained  at  O.N.E.R.A.,  and  referenced  by 
Margason  and  Fearn  [9],  clearly  illustrates  an  entrainment  mechanism  on  the 
wake  side  of  the  jet,  with  the  indication  of  localized  reversal  of  the  cross- 
flow  velocity  component.  In  the  farfield.well  downstream  of  the  plane  Z  *  0, 
and  referenced  in  the  curvilinear  coordinate  system  X^,  Figure  1,  the  char¬ 
acteristic  "horseshoe"  isovel  distribution  for  Uj  and  secondary  vortex 
structure  for  Ug  and  Uj  is  well  documented.  Figure  4  summarizes  the  data  of 
Kamotanl  and  Greber  [10],  as  a  typical  example. 

Much  less  is  known  regarding  the  Interaction  between  the  .iet  and  the 
cross-flow  in  the  nearfield,  i.e.,  in  the  region  0  <  Z/D  <  2.0.  Fearn  and 
Weston's  [11]  correlation  of  experimental  data  on  the  trajectory  of  the  jet 
centerline,  Fiqure  5,  confirms  the  absence  of  data  on  Z/D  <  10.  Fearn  et.al. 
[12]  recently  report  a  compilation  of  limited  data,  for  circular  cross-section 
VSTOL  jets  at  A  »  4  and  A  *  8,  as  obtained  using  a  Dltot-stattc  pressure  rale 
on  the  region  1.0  <  I/O's ,3.0.  Fiqure  6  summarizes  the  measurement  of  the 
velocity  component  rouqhly  parallel  to  the  Z  axis  at  Z/D  =  1.0  and  Z/D  *  3.0. 
These  data  Indicate  some  entrainment  from  the  wake  at  Z/D  *  3.0,  but  none  is 
evident  in  the  data  at  Z/D  *  1.0.  Figure  7  contains  measured  velocity  vector 
distributions  In  the  X-Y  plane  at  the  same  Z/D  stations.  Neither  wake 
entrainment  nor  secondary  vortex  structure  is  Indicated  at  Z/D  *  1.0. 

Limited  entrainment  Is  evident  at  Z/D  s  3.0,  but  the  data  is  too  sparse  to 
confirm  that  the  secondary  vortex  has  been  initiated. 

The  purpose  of  this  study  was  to  expand  and  refine  the  three-dimensional 
parabolic  Navler-Stokes  alqprlthm,  reported  by  Baker,  et.al.  T13],  and  to 
assess  factors  control  lino  prediction  of  the  near-field  flow  development  of 
a  VSTOL  jet  In  cross-flow  at  A  -8.  The  results  of  this  study  are  reported 
herein. 
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a)  Blockage  due  to 
Solid  Cylinder 


b) Blockage  due  to 
Circular  Cross- 
Section  VSTOL 
Jet,  X*  8. 


Figure  2.  011  Flow  Streaklines  on  the  Injection  Plate, 
from  McMahon  and  Mosher  [8], 


Figure  3.  011  Flow  Streaklines  For  A  Circular  VSTOL  Jet, 
From  Margason  And  Fearn  [9] . 
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b)  Transverse  Plane  Velocity 
Distribution,  Xi/D  *  23 


Figure  4.  Experimental  Velocity  Distributions,  Circular  Cross-Section  Jet 
\  ■  8,  from  Kamotanl  and  Sreber  [10], 
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Figure  5  Similitude  For  Circular  Jet  Centerline  Trajectory, 
From  Fearn  And  Weston  [XI] . 


Figure  6.  Experimental  Measurement  of  VSTOL  Jet  Velocity  Distributions, 
Z/D  ■  1.0  and  Z/D  ■  3.0,  Y/D  *  0.,  from  Fearn  and  Benson  [12] 
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THEORETICAL  MODEL 

Three-Dimensional  Parabolic  Navler-Stokes  Equations 

The  3DPNS  equation  set  Is  a  simplification  of  the  steady,  three-dimensional 
time-averaged  Navler-Stokes  equations.  In  Cartesian  tensor  notation  and 
employing  superscript  tilde  and  bar  to  denote  mass-weighted  and  conventional 
time-averaging,  respectively  [14],  the  conservation  equation  form  for  an 
Isoenergetlc  fluid  Is 


1  ■  iq  K1  ' 0 

!i>  ■  fcr  [paiaj  +  psu  +  ^1=7  -  5u]  * 0 
<>  ■  Ijj  5ajk  +  (ck  I  “J  •  “V  It, 


+  pu  j  +  pe  *  0 


L(pe)  *  |jj  pOje  ♦  cj  pu,Uj  |§j  +  e‘  f 


+  c!  £|  .  0 

£  ic 


In  equations  1-2*  p  Is  (constant)  density,  u.  Is  the  mean  velocity  vector, 

J  m 

p  Is  pressure,  is  the  Kronecker  delta.  The  Stokes  stress  tensor  Is 
defined  In  terms  of  the  Reynolds  number  Re  as 


51j  “  (g1j  ‘  I  <51jikk>/Re 


and  pu^uj  Is  the  Reynolds  stress  tensor.  The  fluid  kinematic  viscosity  Is 
v,  and  E^j  Is  the  mean  flow  strain  rate  tensor 
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9U .  3U4 

e  = _ 1  + 

fc1j  '  3Xj  3xj 

Equations  3-4  are  the  transport  equations  for  turbulent  kinetic  energy 

and  Isotropic  dissipation  function,  as  obtained  using  the  closure  model  of 
Launder,  Reece  and  Rodl  [15]  for  the  pressure-strain  and  triple  correlations, 
and 


The  various  coefficients  are  model  constants,  Hanjalic  and  Launder  [16]. 

The  parabolic  Navler-Stokes  equation  set  is  derived  from  equations  1-4 
assuming  the  ratio  of  extremum  transverse  mean  velocity  component  to  down¬ 
stream  (axial)  component  Is  less  than  unity  and  that:  (1.)  the  downstream 
velocity  component  suffers  no  reversal,  (2.)  diffusive  transport  processes 
In  the  downstream  direction  are  higher-order,  hence  negligible,  and 
(3.)  the  overall  elliptic  character  of  the  parent  three-dimensional  Navler- 
Stokes  equation  Is  enforced  through  construction  of  a  suitable  pressure  field. 
Viewing  Figure  1  for  the  VSTOL  jet  problem,  the  x1  (curvilinear)  coordinate 
defines  the  predonrinant  mean  flow  direction  with  scalar  velocity  component  Qj 
of  order  unity,  l.e.,  0(1).  Hence,  assume  0(u2 )  **  0(6)  "*  0(u 3 )  and  0(6) 

<  0(1).  Then,  the  continuity  equation  1  confirms  that  the  downstream 
variation  In  ux  must  be  of  order  equal  to  appropriate  transverse  plane 
variations  of  U2  and  u3;  hence,  for  *  0(1),  *  0(6-1)*  . 


(6) 


Determination  of  the  relative  order  of  terms  In  the  Navler-Stokes  equations 
1-4  Is  straightforward,  cf.  reference  [17].  The  3DPNS  form  denoted  Lp(*)» 
for  the  ui  momentum  equation  Is 


lP(pQi)  =  hi  (p  3i3i>  +  -ii  +  lrr[PUlU2  - 5l2] 

+  1x7  ■  5l3]  a  0 


(9) 
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Should  Xj  correspond  to  a  curvilinear  coordinate  description, the  derivatives 
expressed  In  equation  9  become  the  covariant  derivative.  The  order  of 
pressure  variation  In  the  transverse  plane  Is  determined  by  the  first  order 
terms  appearing  In  equation  2  written  on  U2  and  2$.  Retaining  the  higher- 
order  convection  and  laminar  diffusion  terms  for  generality,  the  30PNS  form 
for  both  transverse  momentum  equations  Is  [17], 


+ 


3f _ 

3xj 3xi 


pVj 


(10) 


which  defines  the  3DPNS  limited  Index  summation  convention,  1  £  (1,j)  £  3, 
2  _<  i  <  3.  Equation  10  Is  an  elliptic  boundary  value  specification  for 
pressure  distributions  In  the  transverse  plane.  The  pressure  field  that 
satisfies  this  Poisson  equation  is  resolved  Into  the  complementary  and 
particular  solutions. 


pUj)  =  Pc(x1)  +  PpUj)  (n) 

The  complementary  pressure  Is  the  solution  to  the  homogeneous  form  of 
equation  10.  The  particular  pressure  pp  satisfies  equation  10,  subject  to 
homogeneous  Dlrlchlet  boundary  conditions  on  boundary  segments  where  p£  Is 
known.  The  30PNS  form  for  equations  3-4  Is 


LP^k*  ‘  3x,  ^u1k*  *  3x( 


p(Ck  f  u,u,  -  5«u)  f£- 

l 


3u 

* 5  +  ■ 0 


LP(e)  ■  ~  (pB,e)  ♦  §^- 


r  Ji  TPTT  1§- 
Lee  U1U*  3x. 


(12) 


Ce 


pu  1  u 


t  e 
i  7 


3ui 

3x7 


0 


(13) 
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With  determination  of  equation  10  as  governing  first-order  phenomena, 

#  — 

It  Is  necessary  to  retain  the  0(6)  transverse  momentum  equations  to 
permit  determination  of  u^.  Retaining  the  highest  two  orders  of  terms,  the 
corresponding  0(6)  transverse  momentum  equations  are 

LS(52k>  ’  IxT  [5“>ak  + 

*  IjjfaA  +p4w'-5w]-° 

Equation  12  employs  the  additional  limited  index  2  <  k  <  3.  The  middle  two 
terms  in  the  second  bracket  are  0(1),  while  the  remaining  terms  are  all  0(6). 
As  a  consequence  of  using  equation  14,  it  becomes  necessary  to  enforce  the 
first-order  effects  of  the  continuity  equation  1  directly  on  Its  solution. 

The  theoretical  concept  Is  to  enforce  a  measure  of  the  continuity  equation 
(solution)  as  a  differential  constraint  on  solution  of  the  0(6)  transverse 
momentum  equations  14.  This  solution  measure  is  chosen  as  the  harmonic 
function  q>(xA),  the  solution  to  the  Poisson  equation 

lP<*>  -  Iq-  <5  M  *  0  <’5 

The  boundary  conditions  for  $  are  homogeneous  Dirichlet  everywhere  at  farfleld. 
Reynolds  Stress  Closure  Model  for  3DPNS 

A  closure  model  for  the  kinematic  Reynolds  stress  -u^uj  ,  appearing  in 
equations  9-13  Is  required.  A  stress-strain  rate  constitutive  equation  yields 
the  kinematic  form  [17], 

-u1uj  3  -kotij  +  c4  *1j  *  CzC^2^1k^kj  +  (16 

where  E^  is  the  symmetric  mean  flow  strain-rate  tensor  given  in  equation  6. 
Retaining  terms  of  the  first  two  orders  of  significance, the  3DPNS  form  for 
the  Reynolds  stress  tensor  Is 


NADC— 81177-60 


^  e'iSXji.sXj  3XZ) 

+  +  lit' 

9X2  (3xi  dXj  j 

.  C-cJCli^ik  +  lHl! 

2  “PjXz^Xs  9X2; 

+  2iitfiit  +  iii* 

^axj^x!  3xj  j 


-  c 


k2  3U; 
*  e  L3X, 


(17) 


Finite  Element  Solution  Algorithm 

The  30PNS  equation  system  has  been  Identified  for  the  dependent  variable 
set  dj(xj)  =  (q)  ■  (uj,  U2,  U3,  p»k,e}T.  The  governing  system  Includes 
equations  1,  9,  10,  12,  13  and  15.  Equations  9,  12  and  13  contain  the 
Initial  value  term  that  permits  a  space-marching  procedure.  Equation  10 
Is  an  elliptic  boundary  value  problem  with  parametric  Initial -value  dependence. 
The  continuity  equation  1  solution  becomes  recast  and  utilized  as  a  differential 
constraint. 

The  generalized  form  for  the  3DPNS  equation  set  Is, 

lP(V  ■  1x7 (5  a*’j>  *  hz  +  ftj]+  sj  ■ 0  <,8) 

where  f^  and  Sj  are  specified  nonlinear  functions  of  their  arguments  as 
determined  by  the  Index  j.  The  solution  domain  Q  is  defined  as  the  region 
l  >  0  In  Figure  1 . 
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The  finite  element  numerical  solution  algorithm  for  equation  16  defines 
the  approximation  q^x^.Ki),  to  the  (unknown)  exact  solution  q^x^,  xi), 
as 

h  M 

qj(xr  xi)  *  qj(xA,  xi)  =  qj(x£,  Xi)  (19) 

Planes  Z  «  C  are  discretized  Into  a  finite  element  distribution  and  in  each 
finite  element  R*,  the  elemental  approximation  Is 

qj*V  Xl^  “  {Nk^xt^T  (2°) 

h 

In  equations  19-20,  j(J)  Is  a  free  index  denoting  members  of  {q  },  and 

sub-  or  superscript  e  denotes  pertaining  to  the  eth  finite  element,  8  =  R* 

t  e  e 

x  xi.  The  elements  of  the  row  matrix  (N^x^)}  are  linear  (k  =  1) 
polynomials  on  x^,  2  <  l  <  3,  spanning  triangles. 

The  finite  element  algorithm  requires  the  generated  error  L^q!))  to  be 

h  ^ 

orthogonal  to  the  functions  {Nk}  employed  to  define  q..  In  addition,  the 

d  h  ^  J 

discrete  approximation  LH(p  )  to  the  continuity  equation  Is  enforced  as  a 
differential  constraint  on  L5(pfl£).  The  theoretical  statement  of  the  finite 
element  solution  algorithm  Is  then 

|R2{Nk}LP(qJ)dx  +  Si  /R2VtNk}LP(ph)dx  =  {0}  (21) 

Equation  21  defines  a  system  of  ordinary  differential  equations  written 
In  the  jet  direction  parallel  to  Xj  in  the  form 

[C]{QJK  +  [U]{QJ}  +  [FLJ]{QL>  +  {SJ}  =  {0}  (22) 

Using  the  trapezoidal  integration  rule, 

AX  i 

{FJ}  =  {QJ}j+1  -  {QJ}j - 2“  +  (QJ}j]  =  (0)  (23) 

and  substituting  equation  22  defines  a  system  of  nonlinear  algebraic  equations 
for  determination  of  the  elements  of  (QJ(xi)}.  The  Newton  Iteration  algorithm 
for  equation  23  Is 
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[J(FJ)]j+1 {5QJ}j*j  -  -  (FJ}P+1 
and 


Additional  details  on  the  algorithm  are  given  In  reference  [17] 
algorithm  Is  operational  In  the  C0M0C:3DPNS  computer  program  [18]. 
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DISCUSSION  AND  RESULTS 


Problem  Specification 

The  analysis  of  a  VSTOL  jet  In  cross-flow,  using  the  3DPNS  computational 
simulation,  requires  specification  of  the  following: 

a)  Solution  domain  extent,  !)  a  RJ  x  xi 

L 

b)  Boundary  conditions  for  q"(x^ )  on  3fl  »  3R  x  xi 

c)  Initial  conditions  for  (x^)  on  *  R2  x  xi  (0) 

d)  Estimation  of  the  x^  pressure  gradient  for  the  u i  momentum  equation. 

Recalling  Figure  1,  the  solution  domain  Is  the  region  Z  >  0,  and  of  sufficient 
lateral  extent  to  permit  specification  of  boundary  conditions  corresponding  to 
the  Imposition  of  cross-flow.  Since  the  plane  Y  *  0  Is  a  plane  of  symmetry, 
only  one-half  of  the  physical  problem  need  be  discretized.  Figure  8a)  Is  a 
sketch  of  the  3DPNS  domain  fl,  denoted  by  vertex  array  A-B-C-D  and  Including 
the  region  Z  >  0.  Since  for  X  *  8,  the  jet  axis  does  not  exhibit  significant 
curvature  on  0  <  Z/D  <  5,  Figure  5,  the  3DPNS  coordinate  system  x^  Is  coin¬ 
cident  with  the  laboratory  reference  frame  X,Y,Z.  The  3DPNS  solution  Is 
marched  parallel  to  the  x3  (Z)  axis. 

Since  the  3DPNS  equation  system  Is  an  elliptic  boundary  wlue  description 
In  the  x2-  x3  plane,  A-B-C-D -In  Figure  8a),  a  boundary  condition  statement 
Is  required  for  every  member  of  qj  (Xj,  x^).  On  D-A,the  symmetry  plane, 
each  dependent  variable  Is  assumed  to  exhibit  a  vanishing  normal  derivative, 
la.,  3q!?/3x2  *  0  on  D-A.  Boundary  segment  A-B  corresponds  to  the  upstream 
cross-flow  boundary,  whereupon  u2  ■  U^-j),  and  k  and  e  are  specified  as  a 
function  of  location  within  the  boundary  layer  assumed  attached  to  the  x2-  x3 
plane.  Both  ui  and  u3  are  assumed-  to  exhibit  vanishing  normal  derivatives  on 
A-B.  Boundary  segment  B-C  Is  assumed  sufficiently  remote  from  the  jet  thatthe 
onset  flow  U^  (-j)  proceeds  unaltered  parallel  to  B-C.  All  other  dependent 
variables  are  assumed  to  exhibit  vanishing  normal  derivatives,  specifically 
Including  03,  the  computed  non-zero  level  of  which  will  correspond  to  deflec¬ 
tion  and/or  entrainment  of  the  cross-flow  by  the  jet  interaction.  Boundary 
segment  D-C  Is  assumed  located  sufficiently  downstream  from  the  jet  such 
that  the  Imposed  cross-flow  Is  recovered,  le  .,  U2(xi,  x^  ■  B-C)  *  U^  (-j). 
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A 


b)  Transverse  Solution  Domain  Discretization  UR*. 


D 


Figure  8.  3DPNS  Computational  Simulation  Solution  Domain. 
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All  other  variables  are  assumed  to  exhibit  vanishing  normal  derivatives 
on  segment  C-D.  Figure  8b)  Is  a  graph  of  the  discretization  R2  *  UR2, 
with  the  cross-flow  distribution  imposed  on  the  boundary. 

The  Initial  condition  specification  for  q^XitOj.x^)  on  ftQ  =  R2  x  xi(0) 
defines  the  basic  problem  and  consists  principally  of  the  jet  cross-flow 
velocity  ratio  and  the  character  of  the  boundary  layer  region  adjacent  to 
the  plane  of  Injection.  As  discussed  at  length  In  reference  [17], It  is 
necessary  to  augment  the  velocity  component  parallel  to  by  a  small 
constant  background  level  to  render  the  3DPNS  description  nonsingular. 

For  a  velocity  ratio  of  X  e  Uj/U^  *  8,  the  background  level  for  ui(xi,  x^) 
=0.2  U.  was  determined  acceptable  and  used  for  these  studies.  The  jet 

J 

Is  defined  to  exist  In  a  semi-circular  region  of  ftQ  centered  on  the  Xj 
axis.  Figure  8.  The  Imposed  cross-flow  u2  Is  assumed  to  exist  only  on  the 
boundary  segments  A-B-C  of  3ft,  starting  at  a  zero  level  at  x.  *  0  and 
proceeding  through  a  boundary- layer  type  evolution  until  u2  »  U^-j)  at 
some  specified  distance  above  the  injection  plane,  say  Xi/D  »  1.0.  The 
Initial  condition  for  u2  everywhere  Interior  to  3ft,  and  for  ua  except  on 
segment  A-B,  Is  self-generated  by  the  3DPNS  analysis  by  the  action  of  the 
continuity  constraint  numerical  algorithm.  The  Initial  distributions  for 
k  and  e  proceed  through  a  similar  boundary  layer  development. 
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Algorithm  Specification 

Recalling  the  ordering  analysis  yielding  the  3DPNS  simplification,  the 
transverse  plane  momentum  equations  were  determined  to  govern  pressure  dis¬ 
tributions  via  a  Poisson  equation.  The  VSTOL  problem  specification  presents 
a  significant  challenge  In  that  the  complementary  pressure  solution  Is  a 
uniform  constant  since  the  Inviscld  flowfleld  exterior  to  the  domain  ft  is 
uniform.  The  particular  pressure  solution  contains  the  entire  pressure 
distribution  resulting  from  the  action  of  generated  Reynolds  stress  and 
transverse  plane  velocity  field  distributions.  However,  no  Initial  speci¬ 
fication  can  be  completed  since  there  Is  no  knowledge  of  the  Reynolds  stress 
and  velocity  field  distributions  on  the  Interior  of  ft  and  near  the  injection 
plane  Z  *  0.  This  is  not  a  critical  Issue  with  the  transverse  plane  velocity 
solution  u^  since  first-order  effects  are  governed  by  the  continuity  equation. 


x 

v 
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However,  the  first-order  momentum  equation  for  ui  contains  the  Xi-pressure 
gradient  as  a  source  term,  and  an  estimate  Is  required  to  achieve  the  first 
solution. 

The  resolution  of  this  issue  was  accomplished  in  this  study  by  definition 
of  a  "virtual  region"  with  an  auxiliary  pressure  gradient  computation,  and 
an  Interaction  procedure  Involving  multiple  overlapping  30PNS  solutions. 

The  "virtual  region"  computation  Is  the  first  3DPNS  solution  generated  on 
the  span  0  <  ;<i/D  <  1.0,  l.e.,  marching  from  the  Injection  plane  a  distance 
of  up  to  one  jet  diameter  parallel  to  the  jet.  In  this  region  It  is  assumed 
that  the  momentum  imparted  to  the  jet  flow  by  the  Injection  process  is  main¬ 
tained  nominally  constant  In  an  overall  Integrated  sense.  The  measure  of 
the  Xi  momentum  Is  the  Integral 
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The  last  form  of  equation  26  Is  the  sum  of  Integrals  over  the  two  dimensional 
plane  of  the  discretization.  Figure  8b).  An  axial  (xx )  pressure  gradient  Is 
computed,  from  the  change  In  momentum  computed  from  equation  26,  according 
to  the  definition 


4L  =  -Mix)  *  -(nij(xi)  -  m0 )  (2 

357  Sj  Ad 

In  equation  27,  A.  Is  the  cross-sectional  area  occupied  by  the  jet  at  the 
Injection  plane,  l.e.,  the  area  of  the  finite  element  domains  interior  to  the 
semi-circle  defining  the  jet.  Figure  8b). 

The  solution  of  equation  27  Is  a  scalar  estimate  of  the  pressure  gradient 
parallel  to  the  xi  axis  required  to  maintain  the  mean  momentum  of  the  jet. 

The  remainder  of  the  30PNS  algorithm  Is  solved, as  discussed  In  the  previous 
section,  specifically  Including  the  particular  pressure  field  which  Is  written 
onto  an  output  file  at  preselected  stations  on  0  <  Xi/D  <  0.5.  In  addition, 
the  entire  virtual  region  30PNS  solution  dependent  variable  set  is  written 
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onto  the  output  file  at  a  pre-determlned  axial  station,  say  Xi/D  -  0.25, 
to  serve  as  a  new  set  of  Initial  conditions  for  the  subsequent  Interaction 
solution. 

The  first  Interaction  solution  overlaps  the  virtual  region  solution  on 
0.25  <  X,/D  _<  0.5,  and  extends  further  downstream  parallel  to  the  jet  axis 
to  Xi/D  a  1.0.  On  the  overlap  region,  the  axial  pressure  gradient  for  the 
ux  momentum  equation  Is  constructed  by  differencing -the  stored  particular 
pressure  distributions  from  the  virtual  region  solution.  The  last  pressure 
gradient,  computed  In  this  manner.  Is  also  extrapolated  and  used  for  the 
3DPNS  solution  on  the  extended  region  0.5  <_  Xi/D  <1.0.  The  first  Inter¬ 
action  solution  particular  pressure  distributions  are  utilized  directly  In 
the  transverse  momentum  equation  solutions,  as  well  as  written  onto  the 
output  file  at  preselected  axial  stations  for  use  In  the  second  interaction 
solution. 

The  second,  and  subsequent  Interaction  solutions  are  restarted  at 
Xx/0  ■  0.25,  or  further  along  the  jet  axis,  using  the  most  rect 3DPNS 
solution  particular  pressure  distributions  to  rompute  the  pressure  gradient 
for  the  ux  momentum  equation.  Each  subsequent  interacts  solution  Is 
carried  progressively  further  downstream  along  the  jet  axis  than  the  span 
for  which  the  pressure  distributions  have  been  computed  from  the  previous 
solution.  Within  each  solution,  however,  the  current  particular  pressure 
solution  Is  employed  In  a  non-delayed  manner  in  the  transverse  plane  momentum 
equations.  In  this  manner,  Item  d)  of  the  3DPNS  analysis  requirements  Is 
met,  assuming  that  the  pressure  gradient  Iteration  procedure  Is  convergent 
and  yields  flowfleld  predictions  In  agreement  with  experimental  data. 

Circular  Jet  Analysis 

The  principal  requirement  Is  to  verify  the  developed  3DPNS  procedure 
for  prediction  of  the  near-field  flowfleld  evolution  of  a  circular  cross- 
section  VST0L  jet  at  A  a  8.  Figure  9  summarizes  the  virtual  region  solutions 
In  terms  of  computed  velocity  vector  distributions  In  the  half-plane  trans¬ 
verse  to  the  jet  axis  at  Xi/D  »  0.25.  Without  the  imposition  of  a  pressure 
gradient  In  the  Qi  momentum  equation,  the  applied  crossflow  ..nmedlately 
penetrates  Into  the  jet  region.  Figure  9a).  Even  though  the  periphery  of 
the  jet  Is  well  defined,  this  penetration  is  at  variance  with  the  experimental 
data.  Figure  7a). 
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c)  dp/dx  ■  0  and  U2tc)  ■  0 

Figure  9.  Comparison  of  30PNS  Virtual  Region  Solution,  Transverse  Velocity 

Distributions,  Circular  Cross-Section  VSTOL  Jet,  X  *  8.,  Xi/D  *  0.25. 
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Computing  the  X*  direction  pressure  gradient,  according  to  equation  27, 
and  applying  It  yields  the  u&  velocity  distribution  shown  In  Figure  9b). 

There  occurs  a  pronounced  turning  of  the  crossflow  towards  the  jet,  and  a 
stagnation  point  Is  predicted  about  one-half  diameter  downstream  of  the 
jet  boundary  In  the  wake  region.  However,  substantial  u2  velocities  are 
still  predicted  on  the  jet  centerline.  For  comparison,  deleting  the 
pressure  gradient  computation  and  defining  u2  =  0  on  the  centerline  within 
the  jet,  yields  the  prediction  shown  In  Figure  9c), which  differs  In  only  a 
small  way  from  Figure  9b). 

These  results  Indicate  that  the  virtual  region  computation,  with  momentum 
conserving  pressure  gradient,  can  produce  an  Initialization  for  the  subsequent 
3DPNS  solutions  that  Is  acceptable.  The  pressure  gradient  application  Is 
crucial  In  that  without  It  the  crossflow  always  penetrates  the  jet.  A 
sequence  of  additional  computations  confirm  that  the  principal  action  of  the 
pressure  gradient  Is  to  counter-balance  the  strong  turbulent  mixing  on  the 
upstream  face  of  the  jet.  Without  this  balance,  a  large  xi-derlvatlve  of  ui 
results,  1e.,  dui/dxi  «  0.  Since  u3  ■  0  and  3«J3/3x3  a  0  on  the  symmetry 
centerplane  and  within  the  jet,  the  balancing  of  first  order  effects 
controlled  by  the  continuity  equation  must  yield  generation  of  a  significant 
velocity  component  u2.  The  momentum  conserving  pressure  gradient  computation 
serves  to  balance  this  coupling  mechanism  within  the  continuity  equation. 

Figure  10  summarizes  the  3DPNS  virtual  region  u^  solution  on  0.5  £  Xi/D 
<  1.0,  as  generated  using  the  momentum  conserving  pressure  gradient.  For 
this  solution,  the  u2  velocity  component  on  the  symmetry  plane  within  the  jet 
•was  also  set  to  zero.  At  Xi/D  *  0.5,  there  appears  a  nominal  radial  outflow 
within  the  jet  region,  a  wake  stagnation  point  about  0.5D  downstream,  and  a 
beginning  Indication  of  entrainment  from  the  wake.  Figure  10a).  There  is 
little  qualitative  change  In  the  appearance  of  this  transverse  plane  velocity 
distribution  on  0.5  <  Xi/D  <  1.0,  except  for  elimination  of  the  rear  wake 
stagnation  point. 

The  principal  requirement  of  the  solution  of  Figure  10  is  to  provide  a 
particular  pressure  distribution  for  the  Ui  momentum  equation  for  the  first 
Interaction  solution.  Since  the  data  of  Figure  10  exhibit  some  level  of 
qualitative  agreement  with  the  experimental  data  of  Fearn  et  al ,  Figure  7, 
the  first  Interaction  can  be  executed.  For  this  and  subsequent  solutions. 
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a)  Xi/D  *  0.5 


c)  Xj/D  *  1.0 

Figure  10.  Summary  of  3DPNS  Virtual  Region  Solution,  Transverse  Velocity 

Distributions,  u2  =  0  on  Symmetry  Plane,  Circular  Cross-Section 
VSTOL  Jet,  \  -  8. 
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u2  on  the  symmetry  plane  within  the  jet  Is  released  and  computed  assuming 
a  vanishing  normal  derivative  boundary  condition.  Figure  11  summarizes  the 
first  solution  as  transverse  plane  velocity  distributions.  The  boundary  of 
the  jet  Is  clearly  defined  at  Xi/D  *  0.5,  and  the  jet  deflection  action  on 
the  crossflow  is  quite  evident.  Particularly  encouraging  is  the  almost 
purely  transverse  outflow  velocity  u»,  computed  In  the  foreward  half  of  the 
jet  region  on  0.5  £  Xi/D  £  1.0.  This  action  Is  In  quite  good  qualitative 
agreement  with  the  data.  Figure  7a).  Also,  this  solution  gives  almost  no 
evidence  of  a  roll-up  of  a  secondary  vortex  Is  Indicated  by  the  data. 

However,  Intrusion  of  a  relatively  large  component  of  u2  occurs.  Interior 
to  the  jet  at  Xi/D  =  1.0,  which  does  not  agree  with  the  Indications  of  the 
experimental  data. 

Figure  12  summarizes  the  second  Interaction  solution  on  0.75  £  X,/D 
^1.25.  A  rather  substantial  Increase  In  entrainment  action  Is  indicated 
at  Xi/D  »  0.75,  In  comparison  to  Figure  lib),  and  a  nominal  transverse 
outflow  for  Ua  Is  maintained  on  the  jet  interior.  Further,  a  wake  region 
stagnation  point  has  reappeared  which  persists  to  Xi/D  »  1.25.  However,  the 
u2  mass  flux  through  the  jet  core  region  persists  at  Xi/D  *  1.0,  and  Is  even 
more  pronounced  at  Xi/D  *  1.25. 

These  3DPNS  predictions  Indicate  that  the  lack  of  knowledge  regarding 
the  detailed  distribution  of  axial  pressure  gradient  In  the  Ui  momentum 
equation  Is  a  constraint.  Even  with  this  limitation,  however,  a  substantial 
quantitative  assessment  of  a  plausible  nearfield  evolution  of  the  crossflow 
has  been  generated.  In  particular.  Figure  13  Is  a  graph  of  pressure 
coefficient  distribution  Cp  *  1  -  pu|/pwl£  on  the  plane  X3/0  =>  0.75. 

Assuming  that  exterior  to  the  jet  region  this  pressure  distribution  is 
nominally  Imposed  through  the  crossflow  boundary  layer,  the  general  distribu¬ 
tion  of  Isobars  Is  In  qualitative  agreement  with  data  [10]. 

In  particular,  the  extremum  level  on  the  lateral  downstream  side  of 
the  jet  Is  due  directly  to  the  strong  entrainment  velocity  field  predicted 
by  the  30PNS  solution  and  confirmed  by  many  experiments.  Further,  Figure 
13b)  Is  a  plot  of  isovels  of  the  jet  direction  velocity  distribution  iii. 

The  Initially  semi-circular  contours  have  been  noticeably  bowed  laterally  in 
the  direction  perpendicular  to  the  Imposed  crossflow.  This  is  also  In  good 
qualitative  agreement  with  data,  cf.  Figure  4a),  and  Is  solely  the  result  of 
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the  Induced  transverse  plane  velocity  field  u^. 

As  discussed,  the  pressure  field  computed  In  the  3DPNS  algorithm  responds 
to  first-order  to  the  Reynolds  stress  distributions,  coupled  to  lower  order 
with  the  transverse  velocity  field,  see  equation  11.  The  accuracy  of  the 
3DPNS  solution  depends  critically  therefore  on  the  stress  distributions. 
Figure  14  summarizes  the  second  Interaction  solution  prediction  of  the 
Reynolds  normal  stress  distributions  at  Xx/D  ■  0.75.  All  three  components 
exhibit  profile  contours  geometrically  similar  to  the  Ox  Isovel  distribution, 
with  the  peak  levels  coincident  with  the  extremum  gradients  In  fix.  Figure 
15  graphs  the  corresponding  3DPNS  solution  for  Reynolds  shear  stress  distri¬ 
butions  at  Xi/D  *  0.75.  The  extremum  gradients  In  utuT  and  uTuT  are  parallel 
to  the  corresponding  transverse  coordinate  direction  u^,  as  expected.  It 
would  be  highly  Informative  to  have  experimental  data  for  comparison  of 
these  predictions. 


SUMMARY  ANO  CONCLUSIONS 

A  three-dimensional  parabolic  Navler-Stokes  numerical  solution  algorithm 
has  been  analy  ed  for  prediction  of  the  nearfleld  flow  development  of  a  VSTOL 
jet  In  subsonic  crossflow.  The  essential  aspects  of  algorithm  definition, 
with  regards  to  Initial  and  boundary  condition  specifications,  has  been 
summarized.  A  momentum  conserving  pressure  gradient  computation  has  been 
developed  to  complete  the  problem  definition  and  facilitate  problem  Initial¬ 
ization.  A  sequence  of  overlapping  Interaction  solutions  has  been  evaluated 
for  prediction  of  a  circular  jet  at  \  ■  8.0.  The  results  of  the  3DPNS  pre¬ 
dictions  have  compared  qualitatively  with  the  sparse  available  experimental 
data.  It  Is  crucial  that  quality  data  be  acquired  to  permit  quantitative 
assessment  of  the  results  of  this  analysis  procedure. 
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a)  Reynolds  Shear  Stress  XiTuT 
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b)  Reynolds  Shear  Stress  XiTuT 


c)  Reynolds  Shear  Stress  uTuT 

Figure  15.  Second  Interaction  Solution  Reynolds  Shear  Stress 
Distributions.  Circular  Cros$-Sect.1on 
VSTOL  Jtet,  X  1 8,  X i/D  ■  0.75. 
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